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ABSTRACT  -  a  steady-state  one  phase  Stefan  problem  correspond! ng 
to  the  solidification  process  of  an  ingot  of  pure  metal  by  conti¬ 
nuous  casting  with  a  non-linear  lateral  cooling  is  considered  via 
the  weak  formulation  introduced  in  [BKS]  for  the  dam  problem. 

Two  existence  results  are  obtained  for  a  general  non-linear  flux 
and  for  a  maximal  monotone  flux.  Comparison  results  and  the  regu¬ 
larity  of  the  free  boundary  are  discussed.  An  uniqueness  theorem  is 
given  for  the  monotone  case. 
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0.  INTRODUCTIOM 


In  this  paper  we  study  the  one  phase  model  of  the  soli¬ 
dification  of  a  pure  metal  in  continuous  casting  submitted  to  a 
non-linear  lateral  cooling. 

In  the  liquid  phase  we  assume  that  the  metal  is  at  the 
melting  temperature,  which  is  zero  after  a  normalization.  In  the 
solid  phase  the  temperature  0  satisfies  the  heat  equation.  The 
ingot  is  extracted  with  constant  velocity  b,  and  the  liquid  - 
-  solid  interface  (the  free  boundary)  is  unknown  but  steady  with 
respect  to  a  fixed  system  of  coordinates  of  |R  ,  in  which  our 
problem  will  be  studied.  Assuming  that  the  free  boundary  $  is 
representable  by  a  surface  z  =  ii)(x,y),  the  steady  Stefan  condition 
i  s 


(0.1)  ®Z"®X*^X  ■  ~  2  =  'J>(i<,y) 

where  X  is  a  positive  constant  representing  the  heat  of  melting. 

In  the  lateral  boundary  one  specifies  a  non-linear 
flux  condition 

(0.2)  -  a0/3n  =  G(0) 


which  expresses  the  law  of  cooling,  and  may  be  quite  general. 
Namely,  we  shall  consider  a  maximal  monotone  graph  G,  which  may  inclu¬ 
de  a  cooling  process  with  cl imati zati on  as  in  Chapter  1  of 
the  book  of  Duvaut  and  Lions  [DL]  • 

This  model  has  been  considered  in  a  particular  case  by 
Rubinstein  [Ru]  and,  with  a  linear  flux  condition  of  Newton  type, 
by  Briere  [Br]  and  Rodrigues  [R]  ,  via  variational  inequalities 
after  a  transformation  of  Baiocchi's  type.  However  this  approach 
doesn't  work  with  a  non  linear  cooling. 

Since  this  problem  has  some  similarities  with  the  dam 
problem,  we  formulate  it  in  section  1  using  the  method  of  Brezis, 
Kinderlehrer  and  Stampacchia  [BKS].  In  sections  2  and  3  we  prove 


-  3  - 


the  existence  theorems,  first  us i ng  compactness  arguments  and  next 
combining  compacity  and  monotonicity  tecniques  for  the  maximal 
monotone  case. 

In  section  4  we  discuss  comparison  properties  which 
show  that  when  the  extraction  velocity  b  is  small  the  ingot  soli¬ 
difies  immediatly  and  there  is  no  free  boundary.  For  some  type 
of  cooling  and  for  ahigh  enough  velocity  b  one  can  show  the 
existence  of  a  free  boundary.  In  this  case  it  is  shown,  in  section 
5,  that  the  free  boundary  is  an  analytic  surface  and  a  weak  solu¬ 
tion  is  also  a  classic  one,  as  in  the  linear  case  of  [R] . 

To  conclude  this  paper  we  give  an  uniqueness  theorem  for 
the  monotone  case  in  section  6,  using  the  techniqueof  Carrillo- 
-Chipot  [CC]  • 


I.  nATHEflATICAL  FORJiULATION 

3 

Let  denote  a  cylindric  domain  in  jR  ,  in  the  form 
2 

fi=rx]0,H[,  where  r  e  |R  is  a  bounded  domain  with  lipschitz 
boundary  91'  representing  a  section  of  the  ingot  and  H>0  its 
height.  We  denote  r^=rx{i},  for  i=0,H,  the  bottom  and  the  top  of 
the  ingot  respectively,  and  by  =  9rx] 0 ,H  [  i ts  lateral  boundary. 
We  have  9f2=r  UTiUTm'. 

0  I  n 

Considering  z  the  direction  of  extraction,  we  can  for¬ 
mulate  our  problem  in  its  classical  form: 


PROBLEM  (C)  :  Find  a  couple  (0,i|>).  such  that 

(1.1)  0>O  in  n  and  0=0  for  0^z<ij)(x,y)  <  H 

(1.2)  AO=b  0^  for  0^4)( X ,y )<z<H 

(1.3)  0=0  on  Tq  ,  0=h(x,y)>O  on 

(1.4)  -  90/3n  «  g(0)  on 


(1.5) 


(1.5’)  0^  ^  Ab,  if  2=4i(x,y)  =  0. 

In  this  formulation  b  and  A  are  positive  constants, 
h  is  a  given  function,  and  g  will  be  specified  in  the  next  two 
sections.  The  reader  will  note  that  the  condition  (1.5')  is  a 
degeneration  of  the  Stefan  condition  (1.5)  in  the  case  when 
the  free  boundary  4>  can  touch  the  Known  boundary  r^,  where  the  melting 
condition  0=0  is  assumed  by  (1.3). 

Let  us  remark  that  by  the  maximum  principle  it  must  be 
e>0  for  z>4i(x,y).  Denoting  by  x*  the  characteristic  function  of 
the  set  fi^={0>O}  and  integrating  formaly  by  parts,  for  every  re¬ 
gular  function  c,  such  that  c=0  on  and  c>0  on  r^,  from 
Problem  (C)  one  has 


(VG-V^  +  bG  c-Abx'^i;_) 


(vo.yc+bG^c-Abc^) 


-A0+bG^)C  + 


1  I?  ^*^“1 

r,u$ur  <i>ur„ 

I  0  0 


g(0)^+ 


1 


•'  r. 


tt(Ab-0j  + 


ll^(0^<^^+0^(J)y-G^+Ab) 


g(0)c. 


1 


where  t" ^  =  !tt^+(fi^+l  .  Therefore,  following  [BKS]  ,  we  introduce  the 
X  y 

weak  formulation  of  Problem  (C)  : 
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PROBLEM 
such  that. 


1  1 

:  Fin(J  a  couple  (0,x)  e  H  (fi) , 


0^0  a . e .  in 

f2,  0=0  on  Tq 

a  nd  0=  h 

on  r^; 

0<x£l  a.e. 

in  Q  and  x=^ 

where 

0>0; 

(1.8)  (V0*75+b0_C-Abxi;,)  +  g(0)C<O,  for  every 

•'fi  ^  ■’r^ 

C  e  such  that  ^>0  on  and  ^=0  on  • 

If  we  consider  a  more  restrictive cVass  of  test  functions 
one  can  introduce  a  more  general  formulation,  which  we  call 
Problem  (P'),  if  we  replace  (1.8)  by 

(1.9)  (VO-Vj;+b0,5-XbxC_)  +  f  g(0)c*0.  VteH\il):t=0  on  F^Ur^. 

•’fi 

It  is  clear  that  every  solution  of  Problem  (P)  verifies 

(1.9) ,  but  the  Problem  (P')  h^s  more  solutions  than  Problem  (P). 
In  particular,  if 


PROBLEM  (P^)  :  Find  0  verifying  (1.6)  and 


(1.10) 


[  (V0*Vc+b0_5)  +  f  g(0)C«O,  VceH^(f2);  t^O  on  F^UF^ 


has  a  solution  0>O,  by  the  maximum  principle,  one  has  0>O 
in  fi  and  (0,1)  is  a  solution  to  Problem  (P'),  which  may  not 
satisfy  (1.5')(see  Proposition  4). 


2.  EXISTENCE  OF  A  WEAK  SOLUTION 

In  this  section  we  assume  the  lateral  cooling  given  by 
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(2.1) 


-  n  (if)  =  g(if.p(x).0(x)),  X  e 


where  p^O  is  a  given  function  representing  the  cooling  tempe¬ 
rature,  and 

(2.2)  g(X,p,e)  is  a  bounded  Caratheodory  function, 
i.e.,  is  continuous  in  6  e  |R  ,  a.e.  {X,p)  e  x  |R^,  measurable 
in  (X,p)  for  all  6  ,  and  maps  bounded  sets  of  r^xlR^xjR  in  boun¬ 
ded  sets  of  |R. 


Since  the  cooling  process  is  determined  by  p,  we  shall 


assume  that 


(2.3) 


(2.4) 


g(X,p,e)<0,  a.e.  (x,p,e)  e  r^xjR^xjR 
g(X,p,e)=0  for  l6l^P.  a.e.  XcT^ . 


such  that 


Consider  a  parameterized  family  of  functions  Xg.eC  {|R) 


,  for  t<0 


(2.5) 


X,(t)  = 


0  1  x^(t)ii 


,  for  0<t<E 


,  for  t>e 


and  so  it  approaches  the  Heaviside  function  when  c'sO. 

Introduce  now  the  following  penalized  problem,  where  for  the  sake  of 
simplicity  we  denote  g(X,p(X)  ,0(X) )  by  g(0)  : 

PROBLEM  (P^)  Find  0^cH\n)  n  C°(R) ,  such  that. 


(2.6)  0^-0  on  r^,  0^-h  on  r^, 


7 


(2.7)  [V0^.Vc+b0^C-AbX(0^)c  1  +  )  g(0^)i:=O,  VtcH^fi):c=0  on  r  ur. 

■'n  .  ^  0  ” 

Assuming  the  functions  h  and  p  verify 

(2.8)  0  <  h(x,y)<M,  a.e.  (x,y)  e  Fj^, 

(2.9)  0  <  p(X)  <  M,  a.e.  X  e  F^  . 

one  can  prove  the  following  "a  priori"  estimate: 

LEMf lA  1  If  0.^  is  a  solution  to  Problem  (P^)  with  assumptions 
(2.2-4)  and  (2.8-9),  one  has 

(2.10)  0  £  0^(X)  <  M,  for  all  X  e  fi  and  0<e5M. 

Pioo^  :  Let  in  (2.7).  One  has 

0  =  {ve^-v[0^]-  +  b0^[0^]”  -  xb  x^(e^)  [©'']■}  +  g(0^)[e^]' 

Jq  2  E  z  jj, 

<  -  {|V[0^]'|^  +  b[0^]'[0^]‘}  =-[  1V[0^]'1^ 

from  which  it  follows  [0^]"=0  and  0^^O. 

From  (2.4)  (2.9)  and  (2.5),  one  has  respectively 

g(0^)  [G^-Mj'^^O  and  Xj.(0'^)  [0^-M]z  =  for  0<e<M. 

Then  c=[0^-M]'*^  in  (2.7)  implies 

0  =  {70^»V[0^-M]'‘'  +  b  0^[0^-Mj'^  -  Xb[GF-M]y} 

Jn  ^ 

.  f  |V[0^-M]^1^ 

■’fi 
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and  therefore  [0^-M]^=Ov,  The  lemma  is  proved.  | 

We  shall  need  the  l“  and  the  Holder  estimates  due  to 
Stampacchia  [S]  for  the  following  eliptic  problem  with  mixed 
boundary  conditions: 

(2.11)  -  Au+bu^=f  in  Q  ,  -1^=9  on  and  u=h  on  r^Ur^. 


LEHNA  2  [S]  The  unique  solution  of  (2.11)  verifies 

(2.12)  llu  11  „  <  C,(  llfll  +  iigll  +  llhll  „ 

L  (fi)  ‘  w  '’P(fi)  LP(r^)  L  (r^ur^) 


(2.13)  Hull  0  ^  <  c.  (llfll  .1  D  +  llgll  Q  +lihll  0  1  _  -  ) 

c°’“(fi)  2  w  '*P(n)  Lp(r,)  c°’'(r^urH) 


for  all  p>3  and  q>2,  and  for  some  constants  Ci  ,C2>0  and  0<a<l 
which  are  independent  of  f,g,h  and  u. 


Pxcoj  ••  See  the  results  of  §5  of  [S]  or  a  more  explicit  result 
extended  to  variational  inequalities  in  Section  2  of  [MS]  | 


Now  we  can  state  an  existence  result  for  the  penalized 
problem,  from  which  we  shall  construct  a  sequence  of  functions 
converging  to  a  solution  of  Problem  (P). 


PROPOSITION  I  Under  assumptions  of  Lemma  1,  and  if 
(2.14)  h  e 

then  there  exists  a  solution  0*"  to  Problem  (P^)  for  all 
0<e<M  satisfying  the  estimate 


(2.15) 


where  the  constants  C>0  .j,and  0<a<l  are  independent  of  e* 


P'too^  :  For  T  e  Bn={T  c  :  HxH  _  <R},  (R>0), 

^  C°(C)“ 

def i ne 


©  =  S^(t) 

as  the  unique  solution  of  the  following  mixed  linear  problem 


0=0  on  r 


0=h  on 


( V0'Vc+b0  ;)=Xb 


J  r 


g(T);,  VccH.  (f^)  :  c  =  0  on 


Since,  by  definition,  0  Xj.  £  1  and  g  is  bounded 
independently  of  t  (for  |t(X)(>M>^  p(X)  one  has  g(X,p(X)  ,T(X)j=0) 
by  (2.4))  one  can  apply  Stampacchi a ' s  estimate  (2.13).  Therefore, 
there  exists  C  >  0  and  0<a<l,  independent  of  t  and  c  such  that 


_  <  C2(Xb+||g 

(fi)  ^  . 


)  1  c 


and  for  R>C  one  has  S  (Bj,)  cz  Bp. 

From  the  compactness  of  the  imbedding  C° ’°‘( fi)‘^C°(n) 
one  finds  that  5^  is  a  continuous  and  compact  mapping  of  Bj^  into 
itself.  By  the  Schauder  fixed  point  theorem  there  exists  a 
function  0^  c  B„  satisfying  0^  =  S  (0^)  ,  which  is  clearly  a  solution 
to  Problem  ( )  • 

The  estimate  in  H^fi)  is  classical,  since  and  g(e^) 
are  bounded  independently  of  c.  | 


THEOREM  1  Assuming  (2. 2, 3, 4)  and  (2.8,9,14)  there  exists  a 
solution  (0,x)c  [h\ji)OC°’“(17)]x  L*(fi)  to  Problem  (P). 


I 


Pioo^  :  By  (2.15)  one  ap  consider  a  sequence  of  solutions  0^ 
of  Problem  (P^),  such  that,  when  e+0 

(2.16)  ^  0  in  H^(ft)-weak 

(2.17)  G^(X)-*'G(X)  uniformly  in  X=(x,y,2)  c  Ti 

(2.18)  X(-(e^)  — -  X  ill  L“(G)-weak  *, 

where  0  is  some  function  belonging  to  ( n)  C® ’‘^(Ti)  satisfying 
(2.10)  and  0_<xil  •  Moreover  in  the  open  set  {0>O)  one  has 
X^(G^)  -*■  1  and  therefore  x=l  ^.e.  in  {0>0}* 

Let  5  e  H^n),  C2.0  on  and  ;=0  on  r^* 

one  has 

<  0 

and  in  the  limit  we  obtain  (1.^8).  The  proof  is  complete.  | 


By  the  Green's  formula  and  since  30^/3n<O  on  F  , 


[y0^.vc+b0^?-Abx  (0^)t  ]+ 


9{0")c= 


1 


10_ 
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3.  THE  CASE  OF  A  fiAXIflAL  nONOTONE  COOLING 

In  this  section  we  consider  the  existence  of  a  weak 
solution  with  a  lateral  cooling 

(3.1)  -  c  G(G)  on  r.,, 

where  G  denotes  a  maximal  monotone  graph  ,  that  is,  G  is  a 
multivalued  function  which  graph  is  a  continuous  monotone  increa 
sing  curve  in  |R^  (see  [B]  ) .  We  shall  assume 

(3.2)  G(0)c=].«,o] 


(3.3) 


[0,+“>[  c:  Oom  (G)={xe  piG(x)^|2}  • 


- 11 


The  weak  forrauVation  of  the  corresponding  problem  takes 
now  the  following  form: 

PROBLEM  (P)  Find  (0,x«g)£:  ( Ji)  j(L*(  n)xL^  {  r-j )  ,  such  that 

(3.4)  G>0  a.e.  in  SI  ,  0=0  on  r  and  G=h  on 

On 

(3.5)  O^Xll  a.e. in  SI,  x  =  l  G>0; 

(3.6)  I  (V0.V^  +  bG  ;-Abx;_)  +  [  g4<0,  V;eH (SI)  :  c>0  on  F  ,;  =  0  on  1'  ; 

;n  ““0  n 

..  ^1 

(3.7)  g(X)  c  G(e(X))  a.e.  X  e  r.|  . 

We  shall  obtain  a  solution  to  Problem  (P)  as  the  limit 

of  a  sequence  of  solutions  to  Problem  (P)  with  a  non-linear  cooling 
given  by  a  function  g  satisfying  : 

(3.8)  g  is  monotone  increasing,  lipschitz  and  such  that  g(0)<0. 

THEOREfl  2  Assume  (3.8)  and  let  h  e  H  ^^^(  ) ,  h>0. 

Then  Problem  (P)  has  a  solution. 

Ptoo |(  :  The  proof  follows  the  lines  of  the  one  in  theorem  1. 

by  considering  the  penalized  problem  (P  )  with  g  satisfying  (3.8). 

c  2 

The  fixed  point  is  now  constructed  in  L  (SI)  by  means  of  the 
mapping 

l2(S1)3  t  -  5  =  T^(t)  c  V. 

where  V={vch\s2):  v=0  on  F^)  and  ^  is  the  unique  solution 
of  the  following  problem 
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(3.9) 


C  c  V  ,  S  =  h  on  r. 


(vc-vc+b^,c)+  g(s)c  =  Xb  VceV:C=0  on  I 


which  is  a  coercive  and  (strictly)  monotone  problem  in  V  by 
assumption  (3.8)  (see  [L] ) .  Denoting  by  R  some  function  in  V, 
which  trace  on  is  h,  and  letting  ?=?-h  in  (3.9)  one  easily 
finds 

IKil  1  <  C  =  C(h)  , 

H'(n) 

where  C  is  a  constant  independent  of  t  and  e* 

1  2 

Since  the  imbedding  H  (fi)  ^  L  (SI)  is  compact, the 
Schauder  fixed  point  Theorem  assures  the  existence  of  a  solution 
to  Problem  (P^.)-  As  in  Lemma  1  one  finds  that  since  g  is 

monotone  increasing  and  g(0)_<0,  and  therefore  one  has  g(0^)’ 

.[0^]'  <  0. 

The  passage  to  the  limit  as  e40  is  straightforward 
since  ©■'— ^  0  in  H  (n)-weak  and  g  is  a  lipschitz  function.  | 


RE^iARK  1  Since  g  is  lipschitz,  by  Sobolev  imbeddings  one 
has  g(0)  e  H^^^(r.| )  ^  L^(r.|)  (see  [A,  p.  218])  and  therefore 
applying  Lemma  2,  it  follows  that 


if  h  e  L  (r^)»  then  0  e  L  (fi) ;  and 

if  he  C°’\r„),  then  0  c  C°’^(Ii),  for  some  0<a<l  .  | 


Since  G  is  a  maximal  monotone  operator  one  can  intro¬ 
duce  the  Yosida  regularization,  defined  by 


where  J^=(I+6G)"^  is  the  resolvent  of  G.  Consider  T*Jg(0), 
that  is  06(I+6G)(t)  .  From  the  roonotoni ci ty  of  I+6G  and  using 
assumption  (3.2)  one  finds  t^O.  Therefore  g^ ( 0) =- ( 0) /6  ^0, 

which  means  that,  for  each  6>0,  the  Yosida  regularization  g^ 
satisfies  the  condition  (3.8)  (see  [B] ) .  So  we  may  apply  Theorem 
2  to  conclude  the  existence  of  a  solution  (0  ,x  )  e  H  (ri:)xL  (Q) 
to  Problem  (P)  with  lateral  cooling  given  by  g^.  We  shall  obtain 
a  solution  to  Problem  (P)  by  considering  a  subsequence  6+0. 

THEOREfl  5  The  Problem  (P)  with  a  maximal  monotone  graph  G 
satisfying  (3.2)  and  (3.3)  ,  and  with  h  e  H^'^^(r^)n  L°°(r^)  has  a 
solution  (Q.x»9)  ^  [H^  (fi)  n  L“(fi)]xL°°(fi)xL'”(r-j )  . 

Moreover,  if  h  e  C°’\r^)  one  has  G  e  for  some  0<a<l. 

Prop ^  :  Consider  the  (unique)  solution  G°  of  the  following  mixed 
problem. 

f  0°eH\fi),  0°=O  on  F^,  0°=h  on  F^ 

i 

(3.10)  ' 

(V0°'Vc  +  b0°c)  +  f  g°(0)c  =  0,  V;eH\^2),?  =  0  on  F  UF  . 

where  g°( t) =Projg^ ^ j 0  is  the  smallest  (in  norm)  number  of  G(t). 
Since  g°(0)£0  it  is  easy  to  show  that  0°^O.  Since  h  c  L  (F^)  one  has 
0°eL'”(f^)  by  (2.12),  and  we  assume  that  0°5M°  =  M^(  h  ,g°(  0) )  . 

X 

Then,  for  every  solution  G  to  Problem  (P)  with  g^,  we 

have 

(3.11)  0  <  0*^  <  0°  <  M°. 

Indeed  (3.11)  follows  by  a  comparison  argument: 
take  in  (1-8)^  and  in  (3.10);  one  has 


w  - 

(3.12)  f  |v[o‘^-e°J*|^,Xbf  [95(G‘^)-g®(o  )]  [0*^-0°] "<0 

Since  0*^j>0  and  x*^=l  in  {g'^>0},  the  middle  terra  in  (3.12) 
vanishes; us  i  ng  g^(0)<0,  together  with 

(3.13)  \9^{t)\  <  |g°(t)i  (see  [Bj  .  p.2D) 
in  order  to  deduce  the  chain 

g^(G^)  >  g^(G°)  >  g^(0)  >  g°(0). 

one  finds  that  the  last  term  in  (3.12)  is  non-negative,  which 
proves  (3.11). 

Using  again  (3.13),  by  (3.11)  one  has 

(3.14)  195(0'^)!  <  lg°(G‘^)i  <  max  [|g°(0)gg°(M°)|]=  £, 

from  where  we  easily  conclude 

llO^II  ,  £  6(  =const .  i  ndepend  ,  of  6). 

H  (n) 

It  follows  that  there  exists  a  subsequence  6+0 

such  that 

(3.15)  G*^  --  0  in  H^(fi)-weak,  and  0<G<M° 

(3.16)  “  X  in  L“(s7)-weak  *, 

(3.17)  g  in  L*'(r^)-weak  *,  with  |ig||^„  5  i  • 

X 

Since  one  can  also  consider  G  -^0  uniformly  in  each 
compact  subset  Kcfl,  one  has  x=i  in  the  open  set  {0>O)* 

Using  the  compactness  of  the  trace  mapping,  one  can 
consider  g'^-*0  in  L^(r.|  )-strong  and  from  (3.3)  in  L^(r^). 

Since  g^(G'^)e  G(J^(g'^)),  it  follows,  by  a  classical  argument 


( [B] ,p.27) .  that  g  e  G(Q)  . 


-  IS 


If  we  assume  h  e  by  Lemma  2  one  easily 

concludes  that  0  e  for  some  0<a<l.  The  proof  is  complete. 

I 


REf'lARK  2  Assuming  that  there  exists  some  vj>0  such  that 
0  c  G(v),  one  can  find  a  more  simple  estimate  in  L“(n)  for  every 
solution  G  to  Problem  (P)  : 

0  <  M  =  max  (v,  llhll  „  )  . 

L  (r„) 

Indeed,  it  is  sufficient  to  consider  in  (3.6)  and  to 

recall  that  the  monotonicity  of  G  implies 
g>0  if  G>M.  I 

REflARK  3  The  results  of  this  section  can  be  easily 
extended  to  the  case  of  a  lateral  boundary  condition 

-  -|2  (X)  e  G(z,G(X))  ,  for  X=(x,y,2)  e 

where  ,  for  each  z  c  ]0,H[,  G(z  ,•)  denotes  a  maximal  monotone 
graph  satisfying  (3. 2), (3. 3)  and  i  in  (3.14)  being  uniformly  bounded  in  z. 

An  interesting  case  could  be  a  lateral  boundary 
submitted  to  N  differents  cooling  zones,  that  is,  when,  for 
i=l,...,N, 


G(z,*)*G.(-).  0=Zjj< 


.<z 


i-1 


<Z<:z  .  <  . 


I 


i\.  COMPARISON  RESULIS. 


If  the  cooling  is  given  by  a  monotone  function  one 
can  adapt  the  technique  of  [BKS]  to  prove  the 
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PROPOSITION  2  Let  0^-  (resp.  )  a  solution  to 

Problem  (P^)  and  correspond! ng  to 
g  and  h  (resp.g  and  fi),  where  g  and  g  areraonotone  functions 
satisfying(3.8)  .  Then  if  fi>  h  and  g<g  it  follows  that  6''>o‘'. 


Set  f ,  ( t )  =  [l -6/ 1] ,  t  e  1r  and  6>0. 


From  (2.7)  and  denoting  one  has 


VrTVC  =  b  {n  +  A[x.(0^)-X^(S^)]}C,  - 

n  Jfi  ^  ^  ■'r. 


[g(o  )-g(o  )]c 


for  every  c  c  =  0  on  F^urj^.  In  particular,  for  ;  =  f^{n), 

which  is  different  from  zero  if  where  9(0^)>g(0^)^g(0^) » 

it  follows 

(4.1)  I  j  VrrV  f^  (n)  1*1  b  L^j  hi  •  I  (n)]  2  I . 


being  the  Lipschiz  constant  of  t  t+A  Xj.(t). 

As  in  [BKS]  ,  (4.1)  implies,  for  any  6>0, 

[  ilog  (1+  )  I  ^  1  C(  =  const .  i  ndepend .  of  6) 

J  n  ° 

from  which  it  follows  0^-0^=  n  1  0.  J 

REMARK  4.  This  argument  also  proves  the  uniqueness  of  the 
solution  of  the  Problem  (P^)  when  g  is  monotone.  Of  course  if 
O(resp.O)  is  a  solution  of  (P)  which  is  the  limit  of  the  subse  - 

E  '  -E  ' 

quence  0  (resp.0  )  the  above  proposition  implies  that  0>0. 


Next  we  shall  prove  comparison  results  with  respect 


to  the  extraction  velocity  b. 
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PROPOSITION  5  Assume  that  there  exists  constants  y,M  such 
that 

(4.2)  0  <  y  ^  h(x,y)  <  M,  a.a.(x,y)  e  . 
and  that  the  function  g  verifies  (3.8)  with 

(4.3)  {t  :  g(t)  =  0}  c 

or  else  that  g  verifies  (2. 2, 3, 4, 9)  .Then  if  ii£qlog(l+^)  a 
solution  G  to  Problem  (P-j)  is  also  a  solution  to  Problem  (P) 

wi th  x= 1  • 


:  If  g  satisfies  (3.8),  then  the  Problem  (P^)  has  a  unique 

solution  (let  X^-O  in  (3.9)).  Moreover  by  (4.3)  one  has  g(0)<O 
(see  Lemma  1 ) . 

Linder  assumptions  (2. 2, 3, 4, 9)  the  existence  of  0 
may, be  shown  essentialy  as  in  Proposition  1,  being  also  g(G)^0, 
by  hypothesis. 

Consider  now  the  function  G^ ( z ) =y ( e^^- 1 )  ( e^^- 1 ) 

Taking  in  (1-10)  and  since  g(0)£O  in  both  cases,  one 

easily  finds  that  0^©^*  Therefore,  if  follows 


3n 


< 


30 


bH 


on 


Using  the  Green's  formula  with  a  smooth  function  c  such  that 

C>0  on  r  and  c«0  on  Tm.  one  has 
—  0  n 


(V0.V;+b  0  ?  -XtnJ.f  9(0)5.  f  (H  ♦  5  <0 

a  "  *  Jr,  Jr„ 
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for  Xb<ijb(e^^-1 ) ”\  This  means  that,  for  all  bH<log(  1 +n/X) , 
(0,1)  is  also  a  solution  to  Problem  (P).  | 

This  proposition  suggests  that,  for  small  velocities 
b,  the  whole  region  is  occupied  by  solid  metal,  since  if  the 
Problem  (P)  admits  only  one  solution  0,  one  has  0>O  in  for 
0<b<p  log(l+y/X).  Conversely  the  next  proposition  suggests  that 
for  big  velocities  the  free  boundary  exists,  since  we  will  show 
that  the  volume  of  the  set  {0>O}  vanishes  when  bt®  • 

PROPOSITION  4.  Under  assumptions  of  the  Theorem  1  or  Theorem  3 


and  denoting  by  the  Lebesque  measure  of  the  set  =  |0{  X)>0}, 

one  has 

(4.4)  1^2+1  1  . 

where  C  is  a  positive  constant  independent  of  X  and  b. 

Moreover,  for  b  big  enough,  one  has 

PjtooJ^.  Let  c=H-z  in  (1.8)  and  in  (3.6),  One  has 

(4.5)  -  0  +b  G-(H-z)  +  Xbf  X  +  f  g(H-z)  <  0, 

•’n  ^  Jn  Jfi  ■'r, 

where  g=g(G)  and  g  e  G(0),  respectively.  In  the  first  case,  g 
is  a  bounded  function  and  from  O<0£M  (see  Theorem  1  and  Lemma  1), 
we  may  assume  with  independent  of  b  and  X.  In  the 

second  one,  by  (3.17)  and  (5.14)  we  have  |lg|l  ^  <z  and  t  is  also 
indepenoent  of  b  and  X  •  ^ 

Denoting  L*  max  (  t  ,£^)  from  (4.5)  it  follows  that 


since  one  has 
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Xh 


' 

•  n 


h  +  L 


(H-z). 


"z  = 


h  and 


H 


a 


e^(H-2)  = 


0  >0 


•'fi 


Recalling  that  0<x£l  and  x=l  in  one  has 
l^+l  .1  XI  |rl  (M+LH^/2)  /Xb. 
which  completes  the  proof  of  the  proposition  .  | 


Now  we  assume  the  existence  of  d,  0<d<H,  such  that 

(4.6)  g{X,p,e)  =  0  for  0<z<d,  V(X,p,e  )er|XjR^xiR 
or,  for  the  monotone  case  (see  Remark  3), 

(4.7)  G(z,.)h0  for  0<z<d<H  . 


Let  (0,x)  ( nesp . (0,x.g) )  a  solution  to  Problem  (P) 
(resp.  (P))  under  assumptions  of  Theorem  1  with  (4.6)  (resp. 
Theorem  3  with(4.7)).  Then  there  exists  6,0<5<d,  such  that 


(4.8) 


0  (x.y.z)  <  Xb[z-d]'^,  V(x,y,z)  e  IT 


(4.9)  0=X*O  "for  0<z<6, 

for  all  b>M/Xd,  where  M^||0|1  is  a  constant  independent  of 

L 

b  (see  (2.10)  and  (3.15)). 

The  proof  of  this  theorem  uses  the  following  lemma. 


.  -  -..-..j'MiNti 


r 
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LEfIMA  3.  Under  assuropti^ons  of  Theorero  4,  one  has 
(4.10)  I  x(Xbx-02)  1  I  (be+Abx-e^)  <  0 


for  0<6<d  and  Z^  =  {(x,y,z)e  JJj  0<z<6}  . 


Pa.0  0^  :  Let  ;=[6-z]'*^  in  (1.8)  or  in  (3.6)  .  One  has 


[-0  +b0,(6-z)+Abx]  1  0, 
J  7  ^  ^ 


because  (4.6)  or  (4.7)  imply  g[6-zj^=0.  Since 


0(6-z)=  0>O  and  O^x^l 

h .  ^  J7.  ” 


it  follows 


x(xbx”0,)  <  (Abx-0_)  <  (b0+Abx-0  )  f  0.  | 

h  ^  ~  J7 .  ♦~JZ-  ^ 


Proof  of  Theorem  A,  ;  consider  y=y(z)=Ab  [z-6]'^  with  6  fixed 
such  that  0<6jcd-M/Ab.  The  function  c=[0-y]^  vanishes  on  z  =  0 
and  for  z^d.  Therefore  g[0-y]'*^=O  and  from  (1.8)  or  from 
(3.6)  ,  one  has 

V0*V[0-y]‘^  +bf  0  [0-y]'^-Ab  x[0”y]^  1  0 

h  J{i  ^  ■'n  ^ 


f  (  lV0l^-Abx0,)  +  f{V0*y  [G-y]'*’-Ab[0-y]^}+bf  0  [0-v]'^  - 
(nsZ^)n{0>O} 


Adding  the  quantity 
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Xb 


^(Xbx-O^)  -  b  j  Xb[0-u]'' 


n\Z, 


"6  "'-6 
which  is  non-positive  by  Lemma  3,  one  obtains 


{02+02+(0^.Xbx)2}+flv[e-p]^i^b|  {0-y)2[0-pr<o, 

^6 


Since  the  last  term  is  zero,  if  follows  that  0<v 
in  n\Z^  =  {z>6}  and  ©^=0^=0.  0^=Abx  in  2^  =  {2<6 }  *5 i nee  0=0 
for  2=0  and  z=6,  one  has  0=0  for  z^6  and  consequently  also 
X=0  for  z<6  •  I 

5,  REGULARITY  OF  THE  FREE  BOUHDARY 

The  goal  pf  Theorem  4  is  to  provide  sufficient 

conditions  in  order  to  assume  the  global  existence  of  a  free 
boundary.  In  this  case  we  shall  prove  that  the  free  boundary 
is  an  analytic  surface. 

We  begin  with  the  following 


PROPOSITION  5.  A  solution  (0,x)  (resp.  (0,x.g))  to  Problem 
(P)  (resp.(P))  satisfies 

(5.1)  -A0+b0^+Xb  Xj  =  0  in  25'(n)  . 

Xz  i  0 


A 


(5.2) 


in  n 


PfLOpj:  The  equation  (B.-'l)  follows  imroediatly  by  taking 
C  ea5(n)  in  (1.8)  or  in  (3.6) 


Choosing  as  a  test  function  in  (1.8)  or  in 
(3.6)  c=rain  (0,en).  where  e>0  and  n  ca(n),ri^0  one  has 
1=  V0*V  min  (0,cn)+b  0  min  (0,eri)-Xb  [min  (0,en)]  <0 


since  x  =  l  in  {0>O}.  Since  min  (0,en)  =  O  on  the  last  integral 
is  zero  and  it  follows 


1  V0l^+e|v0*Vri+b 


(O^en)  {0>en) 


{enO_+0, 
9.  ^  ^ 


[min  (0.en)-erij} 


>  e 


V0* Vn+cb 


{0>en} 


0^n 


-  b 


,^e^[en-0]^ 


from  which  one  concludes 


0,n  <  b  I 


Passing  to  the  limit  e^O,  one  obtains 


(V0*Vn  +  bO  n)  £  0,  V  n  ci6(n):n^0 


and  using  (5.1),  one  deduces  (5.2).  | 


From  (5.1)  it  follows  that  the  function  0  is  locally 


Holder  continuous.  Theref^ore  the  set 

(5.3)  =  {X  e  fi|  0(X)>O} 

is  an  open  set.  Since  x  ■>  s  monotonous  increasing  in  the  z-coordi- 
nate  one  can  introduce 

(5.4)  <l>(x,y)  =  inf  {z  ;  0(x,y,z)>O,  (x,y,z)en} 

where  (I)  is  an  upper  semi -conti  nuous  function,  by  the  continuity 
of  0.  Then  we  can  state. 

THEOREM  5.  For  any  solution  of  Problem  (P)  or  (P)  one  has 

(5.5)  ={0>O}s  {Xefl:z>4){  X  ,y ) } 

where  (f  is  an  upper  semi-continuous  function  given  by  (5.4)  § 

COROLLARY  I.  Under  conditions  of  Theorem  4,  for  all  b>M/Xd, 
one  has 

H  >  (t»(x  ,y  )^d-M/Ad  >  0,  for  all  (x,y)E  T,  which,  in 
particular,  assures  the  existence  of  a  free  boundary.  | 

Consider  now  the  function 

(5.6)  u(x,y,z)  =  [  0(x,y,t)  dt,  for  (x,y,z)  e  fi, 

^  0 


which  is  a  Baiocchi  type  transformation  (see  [BC]  for  instance). 
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THEQREIS  6 .  Let  (0,X)  (•^esp.(0,x.g))  be  a  solution  to  Problem 
(P)  (resp.(P))  under  the  assumptions  of  Theorem  4.  Then  the 
function  u  defined  by  (5.6)  satisfies  the  following  variational 
inequality  in 

(5.7)  u^O,  ( -Au+bu^  +  Xb) >0 ,  u • ( -Au+bu^+Xb)  =  0 , 
and  X  is  a  characteristic  function,  being 

(5.8)  X=x(0)=x(u)  a.e.  in  n 

where  x(v)  denotes  the  characteristic  function  of  the  set 
{v>0}. 


Pfioo  i  :  From  definition  (5.6)  and  recalling  0>0  it  is  obvious 
that  u>0.  Since  0=0^  and  0  satisfies  (5.1)  one  has 

( -Au+bu^+Xbx) 2=  -A0+b0^+Xbx^  =0 

which,  together  with  (4.9)  anc^  O^x^l »  imply 

(5.9)  0=  -Au+bu^+Xbx  <  -Au+bu^+Xb, 

Recalling  (5.5)  it  is  clear  that 

(5.10)  {0>O}  =  {u>0} 


from  which  one  deduces  x=i  i^  u>0,  and  the  third  condition  of 
(5.7)  follows  by  (5.9) . 

From  the  classical  regularity  to  solutions  of  variatio 
nal  inequalities  one  has 

(5.11)  u  e  W^^”(fi)  (see  [KS] ,  for  instance)  and  (5.8) 

I 


follows  easily  from  (5.9)  and  (5.10). 


If  one,  considers  a  linear  fluji 
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(5.12)  g(X.p(X)  ,e(X))  =  a(z)  [e(X)-p(X)] 


with  p>0  and  a(z)=0  for  0<  z <d  and  a(z)=a=const,>0  for 

d<z<H,  then  we  have  that  u  is  the  unique  solution  of  the 
following  variational  inequality  with  mixed  boundary  conditions 
(  see  [Br]  and  [R] ) ; 

u  e  |K  =  {v€H^(I2)|  v^O  in  12,  v  =  0  on 

Vu*V(v-u)+b  u  (v-u)+a  u(v-u)^  h(v-u)-Ab  (v-u)+Q  p(v-u), 
hi  h  J  ■'r^  hi  h^ 

V  V  €  IK, 


where  p( 


fZ 

p(t)dt  for 
d 


z>d . 


In  particular,  this  implies  the  uniqueness  of  the 
solution  of  Problem  (P)  for  a  linear  cooling  given  by  (5.12).  | 

# 

The  transformation  (5.6)  and  its  consequence  (5.8) 
allow  us  to  include  the  study  of  the  free  boundary 


$  =  j2  n 

+ 

in  the  Nnown  results  of  Caffarelli  [c]  Kinderlehrer  and 
Nirenberg  [KN].  In  order  to  apply  these  results  we  must  show 
that  t  has  not  singular  points.  This  may  be  done  by  using  a 
technique  due  to  Alt  [AS,]  for  the  dam  problem. 


imiLJL  Let  Xjje4>  and  B^(X^)c  J2*  Then  there  is  a  cone 
A^c{Xe  |R^|  z<0} .  such  that 

(5.13)  •|^(X)*Vu(X)  .n<0  for  XeB^^2^^o^’  ncA^flS^. 
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Pfiooj :  Recalling  (5. 11). .and  that  U2=0^O  in  n,the  proof  of  this 
lemma  is  a  simple  adaptation  of  Lemma  6.9  of  [KS] ,  page  255,  and 
therefore  we  omit  it.  | 


THEOREfi  7.  Let  (0,x)  { resp . ( 0,x .g) )  be  a  solution  to  Problem 
(P)  (resp.  (P))  under  conditions  of  Theorem  4.  Then  the  free 
boundary  4)  is  an  analytic  surface  given  by 

■S:  2  =4)(x,y)  for  (x,y)  e  T, 

and  0  is  also  a  classical  solution  of  Problem  (C). 


Pxoo{i ;  By  (5.13)  the  function  4)  defined  by  (5.4)  is  a  lipschitz 

function  in  r  and  we  can  apply  Theorem  3  of  [C]  to  conclude  that 
1  2 

(5.14)  <}i  i  s  C  and  ueC  (f2_j^LI4>).  Therefore  from  equation  (5.1) 
and  Green's  formula  one  finds  that  condition  (1.5)  is  verified 
in  every  point  of  the  free  boundary  z=<}>(x,y),  for  all  (x,y)cr, 
by  Corol 1  a ry  1 . 

To  conclude  that  4>  fs  an  analytic  surface  it  is 
sufficient  to  apply  Theorem  1  of  [KN]  ,  using  (5,14)  and  recalling 
that  the  equation  satisfied  by  u  in  has  constant  coefficients.  | 


6.  UNICITY  IN  THE  NONOTONE  CASE 

In  Remark  5  we  have  already  stated  the  uniquenessof  the 
solution  of  Problem  (P)  with  a  particular  linear  cooling. 

Adapting  to  our  problem  the  technique  of  Carrillo  and 
Chipot  ([cc])  we  shall  prove  an  uniqueness  resul  t  for  the  maximal 
monotone  case  assuming  that  x  is  a  characteristic  function,  that 
is,  assuming 

(6.1)  X  =  X(0). 


to  which  we  have  already  stated  sufficient  conditions  in  Theorems 


4  and  6. 


Denote  by  (O^.x^.g^-).  with  and  g.  e  G(0.) 

for  i=l,Z,  two  solutions  of  the  Problem  (P)  and  set 

©o=rain  (G^.e^),  Xo=niin  {xi.X2)*  (<)>1.4>2)- 
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LEfUlA  5 ■  Assuming  (6.1),  one  has 

(6.2)  I  {v(0.-G^) .Vn+b(0.-0^)^n-Ab(x^-Xo)n2  dy  dz 


<  Xb  n(x,y,4)/x,y)  )dx  dy 

JOi 

for  any  n  c  H^fi)  D  C°(fi)  ,  n_^0,  where 


=  {(x.y)  £  r|  (x,y)<(})^(x,y) }  ,  i  =  l,2 


Pa.0  0^  :  Choosing  the  test  functions  +C=+ min(0.j -0^^  ,en) ,  oO,  from 
(3.6)  one  obtains  for  i}<j  (i,j  =  l,2) 


^{V(G.-0.)-V;-.b(G.-G.),C-Xb  (Xi-Xj)^,>-^J^^(gi-gj)c  =  o, 
By  the  monotonicity  of  6,  one  has 


(0i-0o.en)  >  0 


since  it  is  sufficient  to  integrate  in  where 

Then  it  follows 

f 

{V(0^-0jj).V  raTn(0^.-0p,en)+b(0.-0^)2roin(0.-0Q,en) 


-  Xb(Xi--Xo)  [min  (ei-eo.en)]^}  1  0 
or,  using  min  (u ,v)«v- [v-u]'*'. 


28  - 


V(0.-0^).Vn+b 


<  b 


{0.-0  ) 
1  o' 


n- 


0.-0  -,+ 
1  0 


n- 


0..0^n. 


Since  the  x.j  are  characteristic  functions,  integrating 
in  z,  one  has 


■ 

o 

1 

o 

+ 

r  0,-0  ^ 

1  +  r 

r  O-’ 

+ 

(Xi'Xo) 

n-  '  ° 

e  _ 

z 

n-  '  ° 

e 

Vl 

^4 

1 

(x,y,4'.)i 

r\ 

{4).<z<4.^} 


and  (6.2)  follows  by  passing  to  the  limite  c\0  in 


v(0.-0^)  .Vn+bJ  [(o^-Oo)2n-A(x,--Xo)n2]  1 


{0.-0^>en} 


lbJ^(0.-Oo)z 


h - ? — 


+  Ab  ri(x  ,y  ,()>.)  .  | 


THEOREM  8 .  Assuming  (6.1) 
most  one  solution. 


,  the  Problem  (P)  has  at 


P'too |(  :  For  e>0,  consider  a  smooth  function  ct^ ,  such  that, 
O^a^^l ,  and 

a^=l  in  Ap={0^>0}Ur^  and  a^(X)=0  if  d(X,A^)>L. 

Since  l-a^=0  on  {0q>O}  ,  for  all  ncH^n),  one  has 


}  =  0. 


j^{V0^.V(1-a^)n+bOoz(1-aj.)n-XbXp  (l-aj.)n  2 


is  a  test 


For  n  e  H\fi)nc°(n),  n  >  0.  t*(1-a^)n 

function  in  (3.6),  and  it  follows  (since  l-a£=0  on  r^) 


j  {v(0.-e^).v(i-ctjn+b(G.-0^)^ 


(l-a^)n-Xb(x-Xo)  [(l-a^)^]^} 


<  0  (i=l,2). 


Using  (6.2),  we  obtain 

f  <v(0.-Gj.VTi+b(0.-Gj^n-Ab(x-Xo)n,}<  l 


im  xb 
e  +  o 


(“,-n)(x,y,<|)(x,y))=0. 


Choosing  in  this  inequality  n^z  and  n«H-z,  after  a 
simples  calculation  one  obtains 


(x-Xq) 


=  0, 


from  where  one  deduces  X^-~Xq»  'for  i  =  l,2,  which  proves 

the  uniqueness  of  the  solution.  ■ 
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